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PROPULSION  TESTS  IN  LEVEL  ICE  ON  A 
MODEL  OF  A  140-FT  WTGB  ICEBREAKER 

Jean-Claude  Tatinclaux 


INTRODUCTION 

The  United  States  Coast  Guard  initiated  a  mod¬ 
el  experimental  program  in  ice  at  the  U.S.  Army 
Cold  Regions  Research  and  Engineering  Labora¬ 
tory  (CRREL)  on  models  of  the  140-ft  WTGB  at 
scales  of  1:9.273  and  1:24.  The  140-ft  WTGB  is  a 
Great  Lakes  icebreaker  designed  to  operate  in  the 
continuous  mode  of  ice  breaking  at  a  speed  of  l.S 
m/s  (3  knots)  in  0.46  m  (18  in.)  of  level  ice. 

Following  resistance  tests  in  level  ice  on  two 
models  of  the  140-ft  WTGB  icebreaker  (Tatin¬ 
claux  1984),  propulsion  tests  in  level  ice  were  con¬ 
ducted  on  the  larger  model  (scale  1 :9.273).  The  re¬ 
sistance  test  program  was  completed  in  December 
1982,  and  the  propulsion  tests  were  made  in 
April-May  1984.  Between  the  two  test  programs 
at  CRREL,  the  model  was  refurbished  and  was  to 
be  tested  in  ice-free  water  at  the  U.S.  Naval 
Academy,  Annapolis,  Maryland.  This  report 
presents  the  results  of  the  propulsion  tests  in  ice. 


TEST  CONDITIONS 

The  model  tests  were  restricted  to  only  one  ice 
thickness  corresponding  to  a  full-scale  thickness 
of  0.46  m  (18  in.)  and  to  a  narrow  range  of  ice 
flexural  strengths  close  to  a  full-scale  value  of  800 
kPa.  The  majority  of  the  tests  were  conducted  at  a 
model  speed  corresponding  to  a  full-scale  ship 
speed  of  1.54  m/s  (3  knots).  Only  one  test  series 
each  was  conducted  at  model  speeds  correspond¬ 
ing  to  full-scale  values  of  1 .03  and  2.06  m/s  (2  and 
4  knots).  In  addition,  two  resistance  tests  were 
conducted  at  the  equivalent  full-scale  ice  thickness 
of  0.46  m  and  ship  speed  of  3  knots. 


The  model  had  been  refurbished  after  the  resis¬ 
tance  test  program  and  prior  to  the  propulsion 
tests.  The  friction  coefficient  between  the  ice  and 
the  newly  repainted  hull  averaged  0.12,  compared 
to  0. 14  for  the  previous  conditions.  The  lower  fric¬ 
tion  factor  turned  out  to  significantly  reduce  the 
model  resistance  in  ice.  The  results  of  the  friction 
tests  are  given  in  Table  1 . 


EXPERIMENTAL  SET-UP 
AND  PROCEDURES 

The  ship  model  came  equipped  with  a  four- 
bladed  stock  propeller  with  a  nominal  diameter  of 
28  cm  (11  in.)  used  in  earlier,  ice-free  tests  (West 
1975).  The  propeller  shaft  was  connected  to  a 
thrust-torque  dynamometer  rated  at  1 100  N  (250 
lb)  in  thrust  and  11.5  Nm  (100  in. -lb)  in  torque. 
The  input  shaft  of  the  dynamometer  was  driven  by 


Table  1.  Results  of  friction  coefficient  tests. 


N 

(N) 

T 

<N) 

M* 

Top  surface 

783 

108 

0.138 

of  ice 

783 

106 

0.135 

0.134  ±0.005 

783 

101 

0.129 

Bottom  surface 

783 

84.3 

0.108 

of  ice 

783 

86.7 

0.111 

0.108  ±0.003 

783 

82.9 

0.106 

Overall  average:  n,  =  0.121  ±0.014. 

/V  =  total  normal  load  applied  on  ice  sample. 
T  =  measured  tangential  force. 
n,  -  friction  coefficient  =  T/N. 


a.  General  view. 


Motor 


Tow  Post 


c.  Tow  post  assembly. 
Figure  l  (cont’d). 


an  1100-W  (1.5-hp)  variable-speed  motor  via  a 
1:1.7  gear  reducer.  The  shaft  rotational  speed  was 
measured  by  a  magnetic  pickup  over  a  60-tooth 
gear  mounted  on  the  shaft.  The  propulsion  assem¬ 
bly  is  shown  in  Figure  1 . 

In  the  propulsion  tests  the  model  remained  con¬ 
nected  to  the  towing  post  of  the  test  basin  car¬ 
riage.  For  each  ice  sheet  the  carriage,  and  there¬ 
fore  the  model  speed,  was  kept  constant,  and  two 
to  three  runs  at  different  propeller  speeds  were 
made.  The  corresponding  thrust  and  torque  on  the 
propeller  and  pull  on  the  tow  post  were  measured. 
A  positive  pull  indicated  that  the  model  was  re¬ 
strained  by  the  tow  post,  i.e.  that  its  self-propul¬ 
sion  speed  was  higher  than  that  imposed  by  the 
carriage. 

Prior  to  the  actual  propulsion  tests  in  ice,  a  se¬ 
ries  of  forward  bollard  tests  were  made  to  check 
the  overall  propulsion  assembly  by  comparing  the 
measured  pull  against  that  obtained  in  similar  tests 
at  the  Naval  Ship  Research  and  Development  Cen¬ 
ter  (NSRDC)  (West  1975).  These  tests  were  also 
used  to  calibrate  the  dynamometer  against  the 
NSRDC  results,  since  CRREL  does  not  have  the 
equipment  necessary  for  in  situ  calibration.  For¬ 
ward  bollard  tests  were  also  made  immediately  be¬ 


fore  most  of  the  propulsion  tests  in  ice  to  ensure 
that  the  propulsion  and  measurement  systems 
were  functioning  properly. 

In  the  current  tests  the  heave  and  pitch  angle  of 
the  ship  model  were  measured.  The  vertical  dis¬ 
placements  dk  and  d,  at  two  points,  one  in  the 
bow,  the  other  in  the  stern,  at  equal  distances  L 
from  the  center  of  gravity  along  the  ship  center- 
line,  were  measured  by  two  linear  displacement 
transducers.  The  heave  H  and  pitch  angle  </>  are 
given  by 


RESULTS 
Bollard  tests 

The  results  of  the  bollard  tests  are  listed  in 
Table  2.  The  measured  pull  is  plotted  against  pro¬ 
peller  speed  in  Figure  2,  which  also  shows  the  data 
obtained  in  ice-free  water  at  NSRDC.  The  data 
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Propeller 


Test 

no. 

speed  n 
(rpm) 

Pul! 

<N) 

Thrust  7m 
(N) 

Torque  Q,. 
<Nm) 

Thrust 

coefficient  K, 

Torque 

coefficient  10K( 

30 

452 

94.1 

102.9 

3.08 

0.298 

0.319 

558 

145.5 

152.8 

4.81 

0.290 

0.327 

50 

663 

205.8 

219.5 

6.83 

0.277 

0.329 

718 

242.6 

258.8 

8.00 

0.297 

0.328 

772 

280.4 

295.3 

9.27 

0.293 

0.329 

827 

322.4 

339.0 

10.68 

0.293 

0.330 

876 

345.4 

361.7 

11.98 

0.278 

0.330 

772 

279.9 

320.0 

9.26 

0.317 

0.329 

60 

507 

118.2 

134.5 

3.86 

0.309 

0.317 

679 

214.6 

225.9 

7.11 

0.289 

0.326 

843 

332.7 

351.0 

11.06 

0.292 

0.329 

872 

355.9 

375.9 

11.84 

0.292 

0.329 

411* 

791 

2141 

325.2 

9.52 

0.307 

0.322 

421 

837 

326.3 

337.6 

10.96 

0.285 

0.331 

511 

794 

292.8 

297.8 

9.77 

0.279 

0.328 

531 

729 

248.4 

252.1 

8.25 

0.280 

0.328 

611 

857 

342.7 

362.1 

11.49 

0.291 

0.331 

621 

767 

275.0 

286.9 

9.14 

0.288 

0.328 

631 

592 

162.6 

193.9 

5.33 

0.327 

0.322 

•  Tests  411-631  were  run  immediately  before  propulsion  tests. 


Figure  2.  Comparison  between  NSRDC  and  CRREL 
bollard  tests.  Pull  vs  propeller  speed. 
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Figure  3.  Results  of  bollard  tests.  Thrust  and  torque  vs  propeller  speed. 


Table  3.  Results  of  model  propulsion  tests. 


Test 

no. 

tee 

thickness  h. 
(cm) 

Ice 

flexural 
strength  a. 
(kPa) 

Modulus  of 
elasticity 
of  ice  E. 
(MPa) 

Ship 

velocity  V 
(cm/s) 

Propeller 
speed  n 
(rpm) 

Pull 

<N) 

Thrust  Tia 
<N) 

Torque  QIA 
(Nm) 

Heave  H 
(cm) 

Pitch 
angle  <t> 
(deg) 

no 

5.2 

85 

105 

33.5 

790 

7.4 

323.1 

10.97 

0.8 

1.0 

120 

5.2 

85 

105 

33.5 

814 

29.9 

342.5 

11.31 

0.9 

0.9 

130 

5.2 

85 

105 

33.7 

825 

40.6 

352.7 

11.62 

0.8 

0.9 

210 

4.9 

75 

115 

51.7 

870 

18.6 

315.8 

10.81 

0.7 

0.9 

220* 

4  9 

75 

115 

51.7 

0 

-213.1 

— 

— 

0.6 

0.6 

230 

4.9 

75 

115 

52.2 

835 

-8.3 

318.7 

11.40 

0.6 

0.8 

310 

5.3 

76 

112 

66.7 

826 

-91.6 

281.8 

11.88 

0.8 

1.0 

320 

5.3 

76 

112 

66.8 

732 

-137.7 

234.5 

11.86 

0.7 

1.0 

330 

5.3 

76 

112 

66.7 

724 

-132.5 

222.5 

12.06 

0.8 

0.9 

410 

5.3 

75 

146 

50.6 

813 

-51.0 

284.3 

12.03 

0.7 

1.0 

420 

5.3 

75 

146 

50.1 

798 

-45.3 

258.1 

12.06 

0.8 

1.0 

430* 

5.3 

75 

146 

50.2 

0 

-252.8 

- 

- 

0.8 

0.8 

510 

4.8 

90 

76 

50.6 

827 

5.0 

258.6 

11.86 

0.7 

0.8 

520 

4.8 

90 

76 

50.2 

811 

26.1 

263.3 

11.64 

0.6 

0.7 

530 

4.8 

90 

76 

50.2 

775 

-14.6 

228.5 

10.62 

0.6 

0.7 

610 

4.2 

65 

80 

50.5 

869 

90.9 

298.7 

11.69 

0.5 

0.7 

620 

4.2 

65 

80 

50.4 

818 

61.2 

264.4 

10.58 

0.4 

0.6 

630 

4.2 

65 

80 

50.3 

647 

-21.3 

171.7 

7.01 

0.5 

0.5 

•  Tests  220  and  430  were  resistance  tests  (idle  propeller). 


from  NSRDC  and  CRREL  are  in  excellent  agree¬ 
ment.  The  torque  and  thrust  values  are  plotted 
versus  propeller  speed  in  Figure  3  for  the  initial 
bollard  test  results  and  for  those  just  before  the 
propulsion  tests  in  ice;  all  those  results  are  in 
agreement. 


Propulsion  and  resistance  tests 

The  conditions  and  results  of  the  tests  in  level 
ice  are  listed  in  Table  3.  Six  series  of  propulsion 
tests  were  made.  Within  each  test  series  the  ice 
properties  and  model  tow  speed  were  constant, 
and  only  the  propeller  speed  was  varied.  Sixteen 
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propulsion  tests  and  two  resistance  tests  were 
done.  The  results  of  these  resistance  tests  were 
critical  to  the  analysis  of  the  propulsion  tests,  and 
more  such  tests  would  have  been  extremely  useful. 


series  of  propulsion  tests  could  be  estimated  to  be 
83%  of  that  predicted  by  eq  3,  that  is 

=  1 .89  +  0.65  Fi 

■) Bh I 


ANALYSIS  OF  TEST  RESULTS 


Resistance  tests 

From  the  previous  resistance  test  results  (Tatin- 
claux  1984),  the  following  regression  equation  for 
the  total  resistance  in  level  ice  R„  had  been  ob¬ 
tained 


~  =  2.28  +  0.784/^ 
•>  Bn; 


+-  C, [0.042  +  0.0063/%,]  (3) 

7  tin', 


where  •> 
B 
h 
h\ 
l 
e 

c„ 

o, 


specific  weight  of  water 
maximum  beam  at  water  level 
ice  thickness 

ship  velocity 

acceleration  due  to  gravity 

o. /(>>/,) 

ice  flexural  strength. 


The  quantity  R„  is  the  ship  resistance  in  ice-free 
water  and  was  estimated  by 


+  C[0.035  +  0.0052A,,]  +~:  .  (6) 

Of  course  the  variation  in  the  hull  friction  coeffi¬ 
cient  may  affect  the  various  components  of  the  ice 
resistance  to  varying  degrees,  but  the  two  resis¬ 
tance  tests  were  insufficient  to  determine  the  effect 
on  each  component. 

Propulsion  tests 

For  the  conditions  of  the  propulsion  tests  listed 
in  Table  3,  the  corresponding  total  resistance  in  ice 
R„  was  calculated  by  eq  6,  and  the  thrust  deduc¬ 
tion  factor  t  was  then  obtained  from 


/ 


R„  +  Pull 
7m 


(7) 


where  Pull  is  the  average  force  restraining  the 
model  and  T,A  is  the  average  measured  thrust  in 
ice.  The  advance  coefficient 

Jv  =  V/nD  (8) 

the  thrust  coefficient 


R„  =  13.84  V:o2 


(4)  Kr  =  T,A/(Qn2D‘) 


(9) 


for  the  ship  model,  and  by 

R  „  =  586  V:  "  for  V  <  5  m/s 

R  .  =  38  V* for  5  <  V  <  7.2  m/s  (5) 

for  the  full-scale  ship.  In  eq  4  and  5,  R„„  is  ex¬ 
pressed  in  newtons  and  V  is  in  meters  per  second. 

When  applied  to  the  conditions  of  the  two  resis¬ 
tance  tests  (tests  220  and  430),  eq  3  predicted  ice 
resistance  values  (R,  =  R,,-R„.)  of  264.2  and 
289.0  N,  respectively.  The  measured  ice  resistanc¬ 
es,  209.5  N  for  test  220  and  249.4  N  for  test  430, 
averaged  83%  of  the  predicted  values.  In  this  test 
series  the  model  hull  ice  friction  coefficient  was 
only  0.121.  or  86%  of  the  friction  factor  measured 
in  the  earlier  series  of  resistance  tests.  Because  of 
ihe  apparent  correlation  between  the  friction  fac¬ 
tor  and  ice  resistance,  and  for  lack  of  additional 
data,  it  was  assumed  that  the  ice  resistance  (total 
resistance  minus  open  water  resistance)  in  this 


and  the  torque  coefficient 

Kq  =  Qia/(q  n2  D’)  (10) 

were  also  calculated,  with 
n  =  propeller  speed 
D  =  propeller  diameter 
q  =  water  density 

Q,a  =  average  measured  torque  in  ice. 

The  results  of  these  calculations  are  given  in  Table 
4. 

In  view  of  the  scatter  in  the  data  and  the  narrow 
range  over  which  7,  was  varied,  it  was  unjustified 
to  attempt  to  find  variation  of  /,  K,  and  A'g  with 
7, .  Oniy  the  average  values  and  standard  devia¬ 
tions  were  calculated:  t  =  0.214 ±0.150,  K,  = 
0.255  ±0.028  and  =  0.0359  ±0.0025. 

Model  self-propulsion  points 

Two  methods  can  be  used  to  determine  the  seif- 
propulsion  points  of  the  model,  i.e.  the  propeller 
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Table  4.  Analysis  of  model  test  results. 


Test 

no. 

Advance 
coefficient  J , 

Thrust 

coefficient  K, 

Torque 

coefficient  I0KV 

Froude 
number  F, 

Cauchy 
numbei  C. 

Total 

resistance 

in  ice  R' 
<N) 

Thrust 
deduction 
factor  t t 

110 

0.091 

0.306 

0.372 

0.469 

166.6 

246.3 

0.219 

129 

0.088 

0.305 

0.361 

0.469 

166.6 

246.3 

0.198 

130 

0.088 

0.306 

0.361 

0.472 

166.6 

246.5 

0.186 

210 

0.128 

0.246 

0.303 

0.746 

156.0 

222.4 

0.237 

220 

— 

— 

— 

0.746 

156.0 

222.4 

— 

230 

0.134 

0.270 

0.345 

0.753 

156.0 

222  J 

0.327 

310 

0.173 

0.244 

0.368 

0.925 

146.6 

261.0 

0.399 

320 

0.196 

0.259 

0.467 

0.925 

146.6 

261.0 

0.474 

330 

0.198 

0.251 

0.486 

0.925 

146.6 

261.0 

0.423 

410 

0.137 

0.254 

0.384 

0.718 

144.3 

244.1 

0.321 

420 

0.135 

0.239 

0.400 

0.696 

144.3 

242.7 

0.234 

430 

- 

- 

- 

0.695 

144.3 

242.7 

— 

510 

0.131 

0.223 

0.367 

0.737 

191.1 

247.4 

0.024 

520 

0.133 

0.236 

0.374 

0.732 

191.1 

247.1 

-0.038 

530 

0.139 

0.225 

0.374 

0.732 

191.1 

247.1 

-0.018 

610 

0.125 

0.234 

0.327 

0.787 

157.0 

166.4 

0.139 

620 

0.132 

0.233 

0.334 

0.785 

157.0 

166.3 

0.140 

630 

0.167 

0.242 

0.354 

0.784 

157.0 

166.3 

0.156 

•  Total  resistance  /?„  calculated  by  eq  6. 
t  Thrust  deduction  factor  t  calculated  by  eq  7. 


speed  and  corresponding  thrust  and  torque  at 
which  the  pull  is  zero. 

Interpolation  method 

Since  the  series  of  tests  for  each  ice  sheet  was 
conducted  under  practically  identical  ice  condi¬ 
tions  and  model  speeds  but  at  two  or  three  propel¬ 
ler  speeds,  the  self-propulsion  point  of  the  model 
can  be  obtained  by  interpolating  or  extrapolating 
the  data  of  each  test  series.  For  each  test  series  the 
measured  pull,  thrust  and  torque  are  plotted 
against  n\  the  square  of  the  propeller  speed.  The 
straight  line  drawn  through  the  pull  data  points  in¬ 
tersects  the  horizontal  axis  (pull  =  0)  at  the  self¬ 
propulsion  value  of  n\  and  the  corresponding 
points  on  the  thrust  and  torque  lines  give  the 
thrust  and  torque  at  the  self-propulsion  point.  Ex¬ 
amples  of  this  graphical  method  are  given  in  Fig¬ 
ure  4,  and  the  results  of  this  method  of  data  analy¬ 
sis  are  given  in  Table  5. 

Averaging  method 

The  self-propulsion  points  can  also-be  calculat¬ 
ed  from  the  average  values  of  t,  K,  and  KQ  ob¬ 
tained  from  all  the  tests.  For  the  ice  properties  and 
model  speed  of  each  test  series,  the  total  resistance 
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Figure  4.  Graphical  determination  of  self¬ 
propulsion  points  from  test  results. 
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Table  5.  Model  self-propulsion  points. 


Test  conditions 

Interpolation  method 

A  verage  method  * 

Test 

h, 

a, 

V 

n 

T,* 

Cm 

n 

7m 

<2m 

series 

(cm) 

(kPa) 

(cm/s) 

(rpm) 

(N) 

(Nm) 

(rpm) 

(N) 

(Nm) 

100 

5.2 

85 

33.5 

782 

317 

10.8 

852 

313 

12.3 

200 

4.9 

75 

52.0 

846 

317 

11.0 

810 

283 

11.1 

300 

5.3 

76 

66.7 

984 

413 

— 

876 

331 

13.0 

400 

5.3 

75 

51.0 

— 

— 

— 

847 

310 

12.2 

500 

4.8 

90 

50.3 

814 

249 

11.5 

854 

314 

12.4 

600 

4.2 

65 

50.4 

695 

194 

7.9 

701 

212 

8.3 

•  t  =  0.214;  K,  =  0.255;  K0  =  0.0359. 


is  calculated  by  eq  6.  The  required  thrust  is  then 
given  by 

T,*  =  (11) 

1  -/ 

The  corresponding  propeller  speed  is  obtained 
from 


and  the  corresponding  torque  is  given  by 

Q,a  =  Qn1  D'Kq.  (13) 

The  results  of  these  calculations  for  the  six  test 
conditions  are  also  presented  in  Table  5. 

The  results  obtained  by  the  two  methods  show 
significant  differences,  which  is  to  be  expected 
since  even  relatively  small  variations  in  either  t,  KT 
or  Kq  will  affect  the  results.  For  example,  for  test 
series  100  at  a  thrust  of  313  N,  if  one  assumes  a 
value  of  K,  of  0.306  instead  of  0.255  (an  increase 
of  20%),  the  computed  propeller  speed  becomes 
777  rpm  (a  9%  decrease),  while  the  computed 
torque  decreases  by  17%  to  10.3  Nm  and  the  re¬ 
sulting  shaft  power  (PDI  =  2  it  n  Q,a)  decreases  by 
24%  from  1097  to  834  W. 


COMPARISON  WITH  FULL-SCALE  DATA 

Full-scale  trials  of  the  140-ft  Great  Lakes  ice¬ 
breaker  USCG  Katmai  Bay  were  conducted  be¬ 
tween  January  and  March  1979.  The  results  have 
been  reported  by  Vance  (1980a,  b).  The  full-scale 
trial  conditions  and  measured  quantities  of  thrust, 
torque  and  propeller  speed  for  those  runs  made  in 
level  ice  with  no  bubblers  on  are  listed  in  Table  6, 


together  with  the  corresponding  values  of  Jv,  Kr 
and  Kq. 

Most  of  the  full-scale  values  of  the  Froude  num¬ 
ber  F„,  which  ranged  from  0.27  to  2.9,  are  outside 
the  range  of  the  model  Froude  number,  namely 
0.47-0.93.  Similarly  the  values  of  Jv  ranged  from 
0.09  to  0.5  in  the  full-scale  trials  but  only  from 
0.09  to  0.2  in  the  model  tests.  On  the  other  hand 
the  large  majority  of  the  full-scale  values  of  C„, 
ranging  from  155  to  227,  fall  within  the  range  of 
the  model  test  values,  147-191. 

Propeller  characteristics 

The  full-scale  and  model-scale  values  of  KT  and 
Kq  are  plotted  against  Jv  in  Figure  5.  The  full-scale 
thrust  coefficient  decreases  from  about  0.32  at  Jv 
=  0. 1  to  0.24  at  Jv  =  0.4.  In  the  model  test  range 
of  Jv,  the  full-scale  values  of  KT  are  consistently 
higher  than  the  model  data.  The  full-scale  torque 
coefficient  KQ  is  nearly  constant  with  Jv,  with  an 
average  value  of  0.0340,  slightly  lower  than  the 
average  value  of  K0  obtained  from  model  tests. 

The  difference  in  the  model  and  full-scale  values 
of  Kt  can  be  attributed  to  two  factors:  1)  the 
model  propeller  was  a  stock  propeller,  which  may 
have  a  lower  thrust  efficiency  than  the  actual  pro¬ 
peller  used  on  the  USCG  Katmai  Bay;  and  2)  it  has 
been  reported  (Keinonen  1983,  Arctec  Canada 
Ltd.  and  Canadian  Marine  Drilling  Ltd.  1984) 
that  model  ice  floes  broken  in  the  bow  region  by 
model  icebreakers  are  relatively  larger  than  the 
full-scale  floes  and  are  more  readily  entrained  into 
the  propeller  disk,  thereby  reducing  the  thrust  and 
increasing  the  torque  at  comparable  propeller 
speed,  or  requiring  higher  model  propeller  speed 
to  achieve  a  model  thrust  comparable  to  the  full- 
scale  value,  in  addition,  since  the  model  scale  val¬ 
ue  of  Kq  is  slightly  higher  than  its  full-scale  value, 
a  higher  propeller  speed  at  a  given  thrust  will  re¬ 
sult  in  an  excessive  predicted  torque  Qm  and  an 


Table  6.  Full-scale  data  in  level  ice. 


Ice 

thickness  h 
Run  (cm) 

Ship 

velocity  V 
(m/s) 

tee 

flexural 
strength  a ,* 
(kPa) 

Propeller 
speed  n 
(rpm) 

Thrust 

(kN) 

30  Jan  79 

1000 

30.5 

0.72 

610 

148 

86 

1010 

38.1 

2.88 

580 

209 

149 

1020 

27.9 

4.01 

600 

245 

185 

1030 

27.9 

4.52 

580 

259 

196 

31  Jan  79 

1100 

35.6 

0.50 

660 

133 

79 

1110 

36.8 

2.78 

660 

214 

154 

1120 

35.6 

5.44 

650 

252 

191 

1130 

38.1 

4.62 

650 

267 

203 

1200 

38.1 

2.82 

(620) 

204 

145 

1210 

33.0 

3.74 

(620) 

237 

169 

1220 

27.9 

4.51 

620 

263 

196 

9  Feb  79 

1300 

38.1 

1.00 

(640) 

192 

96 

1310 

40.6 

1.87 

(640) 

219 

184 

1320 

41.9 

2.58 

(640) 

247 

213 

1330 

40.6 

2.24 

(640) 

250 

225 

1331 

36.8 

2.87 

640 

254 

218 

m  Torque  Q,„  Advance  Thrust  Torque 

_ (kNm)  coefficient  J,  coefficient  K,  coefficient  IOKv 


24.1 

0.113 

0.312 

0.339 

48.5 

0.319 

0.273 

0.343 

65.3 

0.379 

0.246 

0.336 

75.0 

0.405 

0.234 

0.346 

19.8 

0.088 

0.357 

0.347 

49.6 

0.301 

0.269 

0.334 

67.2 

0.499 

0.240 

0.325 

71.9 

0.401 

0.228 

0.311 

47.5 

0.321 

0.278 

0.353 

60.0 

0.366 

0.242 

0.334 

72.0 

0.397 

0.226 

0.321 

41.9 

0.120 

0.207 

0.349 

53.6 

0.197 

0.305 

0.343 

67.9 

0.242 

0.279 

0.344 

74.9 

0.207 

0.286 

0.368 

75.0 

0.261 

0.269 

0.353 

•  The  values  for  a,  within  parentheses  are  values  assumed  equal  to  the  strength  actually  measured  for  one  run  of  the  corresponding  test  series. 


Figure  5.  Propeller  characteristics.  Kr  and  K0  vs  Jv. 


even  larger  prediction  of  the  required  shaft  power 
Pm,  since  Pm  =  2  r  n  Q,A.  Indeed,  underwater 
video  observation  of  the  model  stern  area  showed 
severe  ice  entrainment  into  the  propeller,  which 
continuously  milled  ice. 

Thrust 

For  the  full-scale  conditions  of  ice  thickness,  ice 
strength  and  ship  speed,  the  total  ice  resistance 


was  calculated  according  to  eq  6,  and  the  required 
thrust  was  obtained  assuming  a  thrust  deduction 
factor  equal  to  that  determined  from  the  model 
tests  (f  =  0.214).  The  predicted  thrust  is  listed  in 
Table  7  «nd  plotted  against  the  measured  thrust  on 
Figure  6a.  It  can  be  seen  that  81%  of  the  calculat¬ 
ed  thrust  values  (13  out  of  16)  lie  between  90% 
and  150%  of  the  measured  values.  The  ratio  of 
calculated  thrust  to  measured  thrust  has  an  aver- 
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Table  7.  Predicted  full-scale  performance. 


Run 

r„ 

<kN) 

From  model  tests  K 

Ku 

From  full-scale  K, 

Ku 

n 

frpm ) 

Cm 

(kNm) 

P,„ 

(HP) 

n 

frpm) 

(kNm) 

Pm 

(HP) 

1000 

112 

188 

40.9 

1080 

169 

31.3 

743 

1010 

194 

247 

70.7 

2452 

236 

61.5 

2038 

1020 

177 

236 

64.7 

2144 

242 

64.4 

2189 

1030 

192 

245 

69.9 

2405 

257 

72.9 

2631 

1100 

144 

213 

52.6 

1573 

190 

39.6 

1057 

1110 

201 

251 

73.1 

2577 

238 

62.6 

2092 

1120 

307 

310 

111.8 

4867 

321 

113.6 

5121 

1130 

279 

296 

101.9 

4236 

297 

97.5 

4067 

1200 

201 

251 

73.2 

2580 

240 

63.2 

2130 

1210 

202 

252 

73.8 

2612 

251 

69.2 

2439 

1220 

199 

250 

72.5 

2545 

261 

75.2 

2756 

1330 

160 

224 

58.6 

1843 

202 

45.2 

1282 

1310 

192 

245 

70.1 

2412 

226 

56.4 

1790 

1320 

219 

262 

80.1 

2947 

246 

67.0 

2315 

1330 

202 

252 

73.6 

2605 

235 

61.0 

2013 

1331 

199 

250 

72.7 

2552 

238 

62.6 

2092 

Table  8.  Average  ratio  between  predicted  and  measured  full-scale  perfor¬ 
mance  (excluding  runs  1120  and  1130). 


Thrust * 

Propeller 

speed 

Torque 

Shaft 

horsepower 

From  model  K,,  K0 

Mean 

1.201 

1.123 

1.345 

1.582 

Std.  dev. 

0.283 

0.175 

0.450 

0.869 

From  full-scale  K,,  Kg 

Mean 

1.201 

1.070 

1.149 

1.263 

Std.  dev. 

0.283 

0.127 

0.293 

0.523 

•  Thrust  calculated  from  predicted  resistance  (eq  6)  and  model  thrust  deduction  factor 
0.214. 


age  value  of  1.20,  with  a  standard  deviation  of 
0.28  (Table  8). 

Propeller  speed,  torque 
and  shaft  power 

Once  the  thrust  has  been  calculated,  the  corres¬ 
ponding  propeller  speed,  torque  and  power  are 
obtained  from  the  thrust  coefficient  XV  and  torque 
coefficient  Kg.  Two  sets  of  calculations  were 
done.  First  the  values  of  Kr  =  0.255  and  KQ  = 
0.0359  obtained  from  model  tests  were  used.  In 
the  second  the  full-scale  values  of  KQ  =  0.0340 
and  Kt  obtained  from  Figure  5  were  used.  The  re¬ 
sults  of  these  calculations  are  listed  in  Table  7  and 
plotted  in  Figure  6b  for  propeller  speed,  Figure  6c 
for  torque  and  Figure  6d  for  delivered  power. 

The  predicted  values  are  usually  higher  than  the 
measured  ones  and  are  even  more  so  when  the  val¬ 


ues  of  Kt  and  KQ  obtained  from  the  model  test  re¬ 
sults  were  used.  The  average  and  standard  devia¬ 
tion  of  the  ratios  of  predicted  value  over  measured 
value  for  propeller  speed,  torque  and  delivered 
power  are  listed  in  Table  8. 


PREDICTION  OF  SHIP 
ICEBREAKING  CAPABILITY 

The  140-ft  WTGB  icebreaker  is  equipped  with  a 
shaft  horsepower  of  2500  hp.  Once  the  equation 
for  the  ship  resistance  in  ice  and  the  propeller 
characteristics  t,  Kr  and  Kg  are  known,  it  is  possi¬ 
ble  to  predict  the  maximum  ice  thickness  h,m.,  that 
the  ship  is  capable  of  breaking  at  full  power  and  at 
a  given  speed.  Such  predictions  were  calculated  us¬ 
ing  the  resistance  equation,  the  value  of  t  obtained 
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Figure  6.  Comparison  between  measured  and  predicted  full-scale  performance. 


from  the  model  tests,  and  the  values  of  Kr  and  KQ 
obtained  from  both  model  tests  and  full-scale 
trials.  The  calculations  were  also  made  for  two  ice 
strengths:  600  and  800  kPa. 

The  results  of  these  calculations  are  presented  in 
Figure  7.  The  model  data  underpredict  the  full- 
scale  performance  for  ship  velocity  below  8  knots 
and  overpredict  it  above  8  knots.  This  cross-over 
occurs  because  the  model  value  of  Kr  is  initially 
smaller  than  the  full-scale  value  for  Jv  <  0.35  but 
is  larger  for  Jv  >  0.35.  Overall,  the  agreement  be¬ 


tween  model  predictions  and  full-scale  data  is  rea¬ 
sonable. 


DISCUSSION  OF  RESULTS 

The  results  presented  in  Figures  6  and  7  and 
Table  8  indicate  that,  on  average,  the  model  tests 
somewhat  underpredict  the  actual  ship  perfor¬ 
mance.  Since  it  is  presumed  that  the  full-scale 
measurements  of  ice  conditions  and  ship  perfor- 


- from  full-scale  values  of  Kr  and  K0 

- from  mode!  values  ofKr  and  K0 


Figure  7.  Predicted  icebreaking  capability  at 
maximum  installed  shaft  power  of  2500  hp. 

mance  are  correct,  the  discrepancy  between  pre¬ 
dicted  and  actual  performance  must  be  due  to  the 
model  test  conditions  and  techniques. 

The  model  test  conditions,  especially  the  Froude 
number,  did  not  cover  the  full-scale  conditions,  so 
the  model  results  could  not  be  directly  extrapolat¬ 
ed  to  the  ship  conditions. 

The  ship  model  was  equipped  with  a  stock  pro¬ 
peller  that  may  nct-have  performed  as  well  as  the 
propeller  installed  on  the  full-scale  ship.  To  check 
this  possibility,  it  would  be  necessary  to  repeat  the 
model  tests  with  an  exact  model  of  the  actual  pro¬ 
peller. 

The  ship  model  hull  may  have  had  a  greater  ice 
friction  factor  than  the  prototype.  The  model 
would  then  yield  exaggerated  predictions  of  ice  re¬ 
sistance  and  consequently  of  the  required  thrust, 
propeller  speed,  torque  and  delivered  power.  Fur¬ 
thermore,  if  the  model  hull  is  rougher,  the  broken 
ice  floes  would  tend  to  “slide”  up  the  hull  less 
readily,  so  more  ice  should  be  ingested  into  the 
propeller,  with  detrimental  effects  on  perfor¬ 
mance.  In  this  respect,  the  full-scale  data  in  Table 
6  indicate  that  the  ship  performance  worsened  be¬ 
tween  the  trials  of  31  January  and  those  of  9  Feb¬ 
ruary  (for  example,  compare  the  results  of  runs 
1331  and  1 1 10).  Also,  the  predicted  thrust,  propel¬ 
ler  speed,  torque  and  power  for  the  later  trial  runs 
are  nearly  equal  to  or  even  less  than  the  measured 
values.  It  is  possible  that  the  hull  roughness  and 


friction  factor  increased  somewhat  during  the  full- 
scale  trials.  The  effect  of  the  friction  factor  can 
and  should  be  investigated  in  model  tests  by  per¬ 
forming  similar  test  programs  at  two  or  more  fric¬ 
tion  factors. 

Broken  floes  in  model  ice  have  been  observed  to 
be  relatively  larger  than  those  at  full  scale.  This 
difference  in  size  has  been  attributed  to  the  two- 
layer  structure  of  the  model  ice  and  to  its  higher 
fracture  toughness  than  that  of  ice  in  the  field. 
Model  ice  may  also  contain  fewer  flaws  and  mi¬ 
crocracks  than  field  ice.  Larger  broken  floes  will 
tend  to  rise  slower  along  the  hull,  because  of  in¬ 
creased  friction,  before  they  reach  the  stern  area. 
More  ice  is  thus  ingested  by  the  propeller,  detri¬ 
mentally  affecting  propeller  characteristics,  name¬ 
ly  the  thrust  deduction  factor,  the  thrust  coeffi¬ 
cient  and  the  torque  coefficient.  The  effect  of  ice 
floe  size  on  propeller  characteristics  and  resulting 
ship  performance  could  be  investigated  by  running 
overload  propulsion  tests  in  a  precut  ice  channel 
for  a  range  of  sizes  of  the  precut  ice  floes. 

The  dimensionless  parameters  F„  and  C„  are  as¬ 
sumed  to  be  the  only  governing  parameters  for  the 
ice  resistance  of  a  given  ship.  The  dimensionless 
ice  resistance,  which  does  not  contain  the  effect  of 
friction  explicity  derived  from  model  tests,  is  then 
assumed  to  apply  to  the  full-scale  conditions.  In 
other  words,  it  is  assumed  that  no  scale  effects  are 
present,  and  that  no  correlation  allowance,  such 
as  is  applied  in  estimating  full-scale  resistance  in 
ice-free  water  from  model  test  results,  is  needed 
for  the  ice  resistance.  Tests  in  ice  on  geometrically 
similar  models  (Geosims)  over  a  wide  range  of 
scaling  ratios  would  be  necessary  to  verify  this  as¬ 
sumption.  However,  tests  in  level  ice  with  large 
models  would  require  very  large  ice  model  tanks, 
which  are  only  beginning  to  be  put  into  operation 
and  are  very  expensive. 


CONCLUSIONS 

On  the  average  the  predictions  based  on  model 
test  results  of  the  ship  performance  compared  rea¬ 
sonably  well  to  that  measured  during  full-scale 
trials.  Among  the  several  possible  sources  of  er¬ 
rors  that  have  been  identified,  duplication  at  the 
model  scale  of  the  ship’s  ice  friction  coefficient  is 
considered  to  be  critical  in  determining  the  ice  re¬ 
sistance  (and  therefore  the  required  thrust  to  be 
delivered  by  the  propeller)  and  the  other  propul¬ 
sion  characteristics,  such  as  propeller  speed, 
torque  and  delivered  power. 


In  spite  of  the  difficulties  in  exactly  modeling 
full-scale  conditions,  the  icebreaking  capabilities 
predicted  from  the  model  tests  are  in  reasonable 
agreement  with  those  observed  during  full-scale 
trials. 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 


Tatinclaux,  Jean-Claude 

Propulsion  tests  in  level  ice  on  a  model  of  a 
140-ft  WTGB  cebreaker  /  by  Jean-Claude  Tatinclaux. 
Hanover,  N.H. :  U.S.  Army  Cold  Regions  Research  and 
Engineering  Laboratory;  Springfield,  Va.:  available 
from  National  Technical  Information  Service,  1985. 
iii,  19  p.,  illus.;  28  cm.  (CRREL  Report  85-4.) 
Bibliography:  p.  13. 

1.  Icebreaker.  2.  Level  ice.  3.  Model  tests. 
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